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1. Introduction 



The AdS/CFT correspondence (or more generally gauge/gravity dualities) provide a recipe 
for computing correlation functions for strongly interacting field theories from classical gravity 
|Q, [|. The analytic properties of field theory Green's functions are constrained by consid- 
erations such as unitarity and causality. For example, causality forces the Fourier transform 
of the retarded Green's function to be analytic in the upper half of the complex frequency 
plane (with our sign conventions). If the AdS/CFT conjecture is true, these analytic prop- 
erties must hold. However, without a proof of the conjecture it is an interesting question to 
ask how these analytic properties arise from the gravitational side of the correspondence. In 
particular, we would like to understand how sum rules for thermal Green's functions emerge 
from gravity. 

Sum rules play an important if sometimes overlooked role in non-perturbative field the- 
ory. There has been recent interest in applying sum rules to investigate the strongly interact- 
ing quark-gluon plasma formed at the Relativistic Heavy Ion Collider (RHIC). To compute 
transport coefficients from lattice QCD calculations, Euclidean Green's functions need to be 
continued to real time. In this difficult game, sum rules provide additional constraints that 
the analytically continued Green's functions must satisfy. To obtain bulk and shear viscosities 
from the stress tensor two-point function, sum rule constraints were employed by for example 
refs. ||, One may also consider sum rules for a conserved current. One famous 

result is the Ferrell-Glover-Tinkham sum rule which relates the "missing area" in the real 
part of the conductivity as a function of frequency to the London penetration depth of a BCS 
superconductor. 1 

In the AdS/CFT context, there have been already a number of papers about thermal 
sum rules. Romatschke and Son [|| derived a pair of sum rules for the stress tensor of J\f = 4 
SU(iV) Super Yang-Mills theory in 3+1 dimensions at large N and strong coupling. Similar 
sum rules for a non-conformal theory dual to a Chamblin-Reall background were constructed 
in [0, |lO|j. The authors of [O] numerically demonstrated an analog of the Ferrell-Glover- 



Tinkham sum rule for a holographic superconductor in 2+1 dimensions, while [12] studied 
an R-charge correlator sum rule for N = 4 Super Yang-Mills theory. In this paper, we will 
rederive the results of || [ll], 12] in a slightly more rigorous and general framework. Our 



point of view is different; we want to see how gravity constraints enforce the sum rules more 
generally rather than checking in specific cases numerically or analytically that the sum rules 
are valid. 

In Section ||, we present the class of gravity theories that we intend to analyze and 
summarize our results. The bosonic correlation functions we study are governed through the 
gravity theory by a single second order differential equation. Instead of considering scalar, 
current, and stress-tensor correlators separately, we unify the analysis through studying the 



1 Recall there exists a Kubo relation between the current-current correlation function and the conductivity. 
The London penetration depth is related to the superfluid density which in turn is proportional to the "missing 
area" . 
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behavior of this single differential equation. In contrast, for the massless fermions, we study 
only correlators of objects with spin 1/2. 

In Section || we review the conditions which lead to sum rules in field theories. Taking 
the corresponding gravity point of view in Sections |] and |5[ we state our assumptions and 
derive sum rules from gravity. In particular, we show that - given our assumptions - the 
Green's functions are holomorphic in the upper half of the complex frequency plane including 
the real axis but not the origin. 

Furthermore, we derive a contraction map for bosons and a contraction-like map for 
massless fermions which we use to bound the large frequency behavior of the corresponding 
spectral functions. From these two major results we derive general sum rules in various 
dimensions and relate them to earlier results. Finally, we discuss our work and suggest some 
future directions in Section ^. 

2. A Class of Holographic Theories & Summary of Results 

Before we describe our class of theories, we present our field theory definition of the retarded 
Green's function for two bosonic operators 2 

G%(t,x)^ie(t)([O\t,x),O\0)}) , (2-1) 

and two fermionic operators 

G R (t,x) = i9(t)({O(t,x),O\0)}) . (2.2) 
Our conventions for the Fourier transform are 

G R (u),k) = [ d d xe^- ihx G R {t,x) , (2.3) 



where we have taken advantage of translation invariance of the underlying theory. We define 
also the spectral density to be the imaginary part of G R (u, k): 

p(u,k) = j.(G R (LJ,k)-G R (Lu,k)^ . (2.4) 

On the gravity side we will not be completely general and instead focus on a particular 
class of d + 1 dimensional space-times with the line element: 

ds 2 = -f(r)e~^dt 2 + — dx 2 + i—. (2.5) 

L 2 f(r) 

We will assume that the space-time contains a horizon at r = where f(rh) = and 
f'( r h) = fh 7^ 0. We also assume that at large r, the space-time becomes asymptotically 
anti-de Sitter with a radius of curvature L: 

f = J? + ... , e-x/ = _ + ... . (2.6) 



2 Note that our definition of the Green's function differs from the one given in |lq] by a minus sign. 
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These space-times are conjectured to be dual to thermal field theories with d space-time 
dimensions and have appeared many times before in the AdS/CFT literature. They are 
sufficiently general to allow us to consider field theories at nonzero temperature and charge 
density and also, if need be, in a superconducting or superfluid phase [|^, 14, 17 1 - 

The prescription for calculating Green's functions was provided in the original papers 
H |3[ and later extended to nonzero temperature by [18, [HJ. The procedure begins by finding 
a solution to a system of differential equations describing fluctuations of classical gravitational 
fields. In this paper, we will consider the two-point functions for a scalar, a conserved current, 



the stress tensor, and a spin 1/2 field in the space-time (2.5). 

For bosonic fields, the following versatile second order differential equation will play a 
big role in this paper: 



V(4>)(r) = (j)"{r) + T^r^'ir) + T o (r)0(r) = , 



where ' denotes d/dr, 



F' n 
T X = — + - 
b r 



u? Y 
and T = - - , 



(2.7) 



(2. 



Y is a potential term specified below, and F = fe~ x ^ 2 . This differential equation can be used 
to study two-point functions of scalar, vector and tensor fields from gravity. 

After a brief review of sum rules in Section ||, in Section |I| we will study the analytic 
structure of the Green's functions for a scalar O, a component of a conserved current J x , and 
a component of the stress-tensor T xy using the differential equation (|2.7j ). We discuss the 
reality and positivity properties of Gu(u,k): Gn(uj,k)* = Gn(—u,—k) and up(uj,k) > 0. 
We show that Gr{oj, k) continued to complex uj, has no poles or branch cuts in the upper half 
of the complex plane including real We discuss the large uj behavior of Gr{uj, 0) from 

which we can justify the sum rules. Finally, we repeat this analysis for massless fermions in 
Section ||. 

Let us first explain how ([T?]) governs the bosonic correlators and describe the results for 
fermions and bosons in greater detail: 



Scalars 

Consider a real scalar operator O dual to a minimally coupled scalar field (j) in the gravity 
dual: 



Sa 



d: 



■d+l. 



-9 [{d^W 



+ m 



2 j.21 



(2- 



For n = d — 1 and Y 



[m 



+ k 2 L 2 /r 2 )e~ x ^ 2 the differential equation ( |2.7| ) is the equation 
of motion of a scalar field (f> of mass m and space-time dependence <fi ~ e -t^t+ikx^ rp^g 
momentum space two-point function for O can be read off from the large r expansion of (j). 
The expansion takes the form 



ar 



A-d 



(l + ...) + 6r- A (l + ...) , 



(2.10) 
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where the ellipsis denotes higher order terms in r and A(A — d) = m 2 L 2 . Note that A 
can be interpreted as the scaling dimension of the dual operator in the field theory. From 
gravity, solving the differential equation with ingoing boundary conditions at the event horizon 
r = rh, up to an overall normalization and possible logarithmic contact terms, we have 
GR(u,k) ~ -b/a. 

We do not press the analysis as far as a demonstration of sum rules, with the exception 
of m = 0, where the analysis overlaps with the stress tensor correlator described below. 

Conserved Currents 

Conserved currents in field theory are dual to gauge fields in the gravity dual. If we set 

n = d — 3, (j) = A x , and Y to a particular form we describe in Appendix [A| (see ( |A.10[ )), then 

a solution to ( |2.7| ) describes the G xx (uj,0) component of the current-current Green's function 

in the dual field theory. Indeed, the presence of Y allows us to study Gr as a function of 

charge density and also follow its behavior through a superfluid phase transition. 

The Green's function can be found from the behavior of the gauge field near the boundary. 

Restricting to d = 3 and 4, the x-component of the gauge field takes the form A x = A x ^ + 
j4 (l) r -(d-2) + ( r -(d-i))_ Being 

more careful about normalization than we were in the scalar 
case, the Green's function is then (see for example fll~8|| ) 

A — 1 ) — 1 Fr d ~ 3 A ; 

where e is a coupling constant for the £7(1) gauge field strength. In d = 4, this expression 
will contain a logarithmically divergent contact term that needs to be regularized. 

In the cases d = 2, 3 and 4 when Y = x = 0, an d / = r 2 (l — (rh/r) d )/L 2 , the equation 
( |2.7D is exactly solvable. This case corresponds to a thermal field theory at zero charge density 
(and also not in a superfluid or superconducting phase) . One can explicitly verify the analytic 



structure of the Green's function. We review these known results in Appendix A.l. 

The more interesting general case where F / we save for Section Given some 
assumptions about the behavior of Gj^(0,0), we are able to prove that the Ferrell-Glover- 
Tinkham sum rule holds from gravity considerations. 

Stress Tensor 

The stress tensor T^ v in field theory is dual to fluctuations of the metric in the gravity dual. 
It is well known that off-diagonal metric fluctuations Sg xy = (f>{r)e~ lult+lkz r 2 / L 2 that are not 
functionally dependent on x or y obey the equation for a minimally coupled scalar of mass 
m = 0. Thus, a solution to ( [2.7D for n = d — 1 allows a calculation of a particular component 
of the stress-tensor two point function: G^' xy (uj,k) where the x and y components of k are 
set to zero. Up to contact terms and regulators we describe in more detail in Appendix [B|, 
the retarded Green's function is 



i — >oo 
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where k is the gravitational coupling constant. This correlator is of particular interest since 
it is related to the shear viscosity transport coefficient [^]. Among other results, we recover 
in general space-time dimension a sum rule proposed by Romatschke and Son ||] for j\f = 4 
Super Yang Mills, i.e. the d = 4 case. For the sum rule, we restrict to Y = and hence 
k = 0. Note that k = is still the interesting case for the shear viscosity computation. 

Fermions 

The operators O, J x , and T xy are all Hermitian. By way of counterpoint, in Section |5], we 
consider sum rules for a non-Hermitian operator — a minimally coupled massless fermion in 
an asymptotically AdS^+i space. The relation G^f(u},k)* = G^(— uj, — k) (with a, (5 spinor 
indices) does not hold. We show p(uj, k) > for real co ^ 0. 3 We prove that all (quasi)normal 
modes must either lie in the lower half of the complex frequency plane, or at u) = 0. Finally, 
we study the large uj behavior of G r (uj, k). These are the ingredients from which we derive a 
sum rule for massless, charged fermions in arbitrary d. For d = 3, the sum rule we find is a 
special case of a sum rule proposed in an Appendix of [^] . 

3. Sum Rules from Field Theory 

Consider a retarded Green's function G r {uj) and the corresponding spectral density, p(uj) 
fl2.4|) . While this section can be generalized to the case where G r (uj) has a multi-index 
structure, for clarity we will assume that G r (uj) is the single component retarded Green's 
function of a scalar bosonic operator. We make two assumptions about the retarded Green's 
function: 

1. G r {uj) is holomorphic in the upper half plane, including the real axis. 

2. limj^i^oQ Gr[uS) = if Imw > 0. 

The first assumption enables us to apply the Cauchy integral theorem: 

„ , , -x / G R {z)dz I G R (z)dz 

G R (u + »e) = f — — ; 0=<f>— — — — . 3.1 

J c lm{z — uj — ie) J c liriyz — uj + ie) 

For w£l and some r G M + , we choose C to be the contour that travels along the real axis 
from — r to r and then along a semicircle in the upper half plane back to — r. Our second 
assumption tells us that the integral along the semicircle must go to zero as r — >• oo. Then 
we have 

n! ... r G R ( Z )dz r G R ( Z )dz 

G r (uj + ie) = km / — — — ; = hm / — — — — . (3.2) 

r->oo J_ r 2tTI{Z — UJ — ie) r->oo J _ r Z7Tl(Z — UJ + l€) 



3 The missing factor of uj in this inequality can be traced to the different statistical distributions obeyed by 
bosons and fermions. The Fourier transform of the variance, (O^ (x)O(0)}, is non-negative, as can be seen by 
decomposing O(x) into a complete basis of eigenstates. The variance is then related to the spectral density 
by a Bose-Einstein or Fermi-Dirac distribution function, depending on whether 0(x) is a bosonic or fermionic 
operator. 
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We subtract the complex conjugate of the second integral from the first integral in (3.2) and 



take the limit e — > to find a "spectral representation" of the retarded Green's function: 

G R (uo) = lim / . (3.3) 

e^0+ J IT Z — UJ — 16 



If assumption 2 is not satisfied, then a modified version of (|3.3[ ) can often be constructed, 

8G R (u) = lim / *_WL. , (3.4) 

e->0+ J IT Z — UJ — IE 

where bp = Im 5G R . The precise definition of 5G R (to) depends on the regularization proce- 
dure. For example, suppose 

G R (u) = 5G R (to) + c + ciw + . . . + c n to n + (d + d-iuj . . . + d s w s ) log(-iw) (3.5) 

where 5G r (uj) does satisfy assumption 2. We take the branch cut of the log to lie along the 



negative imaginary axis. The relation (3.4) then replaces ( |3.3| ). Other regularizations are 
often possible. The coefficients q and dj may be independent of a parameter in the theory, 
for example temperature T. In this case, we could construct a regulated 5G R by considering 
the difference G R (T\) — G R {T%). Another way to regulate the large u divergence is to take 
derivatives: 

d n+l r . , .. fdz (n + 1)1 p(z) 

■Gr(v) = hm / — — — - — — — . (3.6) 



While the entire function Gr(cu) is generically difficult to determine in an arbitrary quantum 
field theory, it is often possible to obtain the q and dj in the limit Imw —> oo, for example 
through an operator product expansion in an asymptotically free theory. 

Given an appropriately regulated version of (|3.3|) , people often call the spectral represen- 
tation evaluated at oj = a sum rule: 

6G R (0) = [ — ^ . (3.7) 

J TT U 

By construction, this object is convergent for large uj, but in discarding the ie regulator, 
we have to be careful about the convergence properties of the integral at uj = 0. From the 
definition of the retarded Green's function, it follows for bosonic Hermitian operators that 
G R (uj,k)* = G R (—L0,—k) (for real to and k). Let us restrict to the case k = 0. In a small 
to expansion of G r (oj, 0), even powers of uj have real coefficients and odd powers have pure 
imaginary coefficients. Provided 5G R (0, 0) is finite, the sum rule should be well defined. 4 
There are also generalized sum rules that involve derivatives of G R (to): 

1 f G *(°) = / ^TT U") ~ t , (3-8) 



(2j)! duj 2 i ny ' J ttuj 2 ^ 1 



=i 



4 We will see a case in the appendix, for the current-current correlation function in d — 2 space-time 
dimensions, where Gr{uj) ~ i/u at small frequencies. In this case, we have to be more creative to write a sum 
rule. Note that this behavior violates assumption 1 above. 
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where 2j > n,s. For these generalized sum rules, the a« are chosen to ensure that the 
integral converges at cj = 0. Note we have taken advantage of the reality properties of a 
small u> expansion of Gr{uj). Some specific examples of this type of generalized sum rule for 
current-current correlators are given in the appendix, ( A.27| ) and ( |A.37 ). 



4. Sum Rules from Gravity 

We would like to see how correlation functions for bosonic Hermitian operators derived from 
the gravity side of the AdS/CFT correspondence satisfy sum rules of the form ( |3.7| ). (The 
fermionic case is postponed to Section g.) Consider a generic retarded Green's function 
Gr(lo, k). We will establish, from gravity, that Gr(oj, k) is holomorphic in the upper half plane 
and also for real We will check that Gr(oj, k)* = Gr{—uj, —k) and that cop(uj, k) > 

(for real uj and k). We will determine the large frequency behavior of Gr(oj, 0) so that we can 



appropriately regularize (|3.3|) . We will not be able to show in general that Gr(0,0) < oo." 



4.1 Imw > 

There is a very good, well known physical reason to expect that Gr(lo) is holomorphic in the 
upper half plane. The reason is field theoretic and related to causality. Consider the Fourier 
transform back to real space 

G R (t) = J ^e-^G R (u) . (4.1) 

If Gr(ui) has poles and branch cuts only in the lower half plane, then the contour can be 
closed in the upper half plane for t < and the integral evaluates to zero, as expected if the 
response of the system is to be causal. 6 

As emphasized in the introduction, we would also like to have a gravitational reason. In 



Appendix D.l we show that there are no branch cuts in the class of theories we analyze here. 
So we consider only meromorphic Gr(uj, k), and the only non-analyticities are poles which ref. 
[ |T^ 1 has argued are dual to quasinormal modes in gravity. The locations wS* of poles of the 
boundary theory Green's function in the complex frequency plane are exactly the quasinormal 
frequencies Wn™ on the gravity side. A quasinormal mode is a solution where the bulk field 
has ingoing boundary conditions at the event horizon of the black hole and where the leading 
behavior of the field vanishes, e.g. = for the current-current correlation function. Given 
the assumed time dependence e~ lujt , a quasinormal mode in the upper half of the complex 
plane is a solution that is exponentially growing in time, indicating an instability of the metric 
and/or matter fields in the gravity dual. An interesting consequence of gauge/gravity duality 



5 See footnote [I| 

6 Note that we could redefine our Green's function such that the integration contour passes above all of the 
poles. While the resulting Green's function will be causal, the contour deformation is non-standard, and the 
poles in the upper half plane indicate instabilities in the field theory. We thank Karl Landsteiner for discussion 
on this point. 
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is thus that the boundary requirement of causality Im(wn° le ) < corresponds to the gravity 
requirement of stability Im(wn nm ) < 0. 7 

Just because a pole in the upper half of the complex plane would be a bad thing doesn't 
mean it can't happen. (Indeed, such poles drive the holographic superconducting phase 
transition described in |TT| , P^D- Given some additional mild technical assumptions, we now 
prove for the class of scalar, current, and stress tensor correlation functions governed by the 
differential equation fl2.7|) introduced above, there are poles only in the lower half of the 



complex plane and at oj = 0. In particular, we make some general assumptions about (2.7). 
We assume that n > 0. Next, we assume F > and Y > on the interval r\ t < r < oo. We 
also assume the near horizon behavior F ~ (r — r^jF^ > and Y^r^) < oo and the large 
r behavior F ~ r 2 and Y ~ r ~ 2A where A > 0. From these assumptions, it follows that 
a quasinormal mode solution to this differential equation has the near boundary behavior 8 



-n-l 



A ~ (r - r h )~ luj/Fh . (4.2) 

The near horizon behavior is chosen to give a retarded rather than an advanced Green's 
function. With an implicit time dependence e~ lujt , the boundary condition corresponds to a 
plane wave traveling into the event horizon. 

Such a differential equation fl2.7|) can be derived from the one dimensional effective action 

(4.3) 



,,2' 

S A = / drFr n 



\A'(r)\ 2 +(^-^)\A(r 



rh 

Assume that there exists a quasinormal mode solution A to (|2.7| ) with frequency oj. Because 
the differential operator is real, there will be a second complex conjugate quasinormal mode 
A with frequency oj. Consider = S A — Sa, evaluated on implicit solutions A and A. We 
integrate by parts, using the equation of motion for A in the first Sa, and that for A in the 
second. We find that 

poo n 

= Fr n [A' A - AA'] |~ + (u 2 - co 2 ) / dr — \A\ 2 . (4.4) 

Because of the large r behavior of A, the first term must vanish evaluated as r — > oo for 
n > 0. Now if oj lies in the upper half of the complex plane, the first term will also vanish 
evaluated at r = rh- In this case, Co 2 = oj 2 and we find that co lies on the positive imaginary 
axis. For uj in the lower half plane, in contrast, the differential operator is not self-adjoint 
and uj 2 need not be real. 



7 Causality can also be studied considering the front velocity, and analyzing the poles of Gi?(oj,fc) in the 
complex momentum plane (for real frequencies) 

8 If A = 0, as for example happens for a massive scalar field, we have Y ~ m 2 near the boundary. Assuming 
m 2 > 0, it follows that A has an even smaller exponent A ~ r -i/2(n+i +v V+i) 2 +4m 2 £ 2 ) at large r The 
argument we are about to present justifying the absence of quasinormal modes in the upper half plane continues 
to hold. 
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Assume now that there exists a quasinormal mode in the upper half plane for which 
u 2 < 0. In this case, Sa is positive definite. If we integrate by parts and use the equation of 
motion fl2.7D , Sa reduces to the boundary term 



S A = Fr n A(r)A'(r)r . (4.5) 

1 ' h 

From the near boundary and near horizon behavior of F and A, Sa = 0. Thus, no quasinormal 
modes in the upper half of the complex plane exist. 

Note if we weaken the condition on Y and allow Y to become negative and develop a 
potential well, then such a quasinormal mode can exist. As an example of such a quasinormal 
mode crossing into the upper half plane, consider the case where A is a charged scalar field in 



the holographic superconductors ]IT], [13, 14|, M. Because the scalar is charged, it does not 
fall in the class of fluctuations which we address in this paper, and its instability is thus not 
ruled out by our argument. As a similar example, consider the holographic superconductor 



involving instability of a neutral scalar [14]. Here a negative mass squared causes Y < 
violating our initial assumption Y > 0, and this case again is not in our class of models. 9 In 
both cases, this scalar instability causes the phase transition to the superconducting state. 
In the new stable phase, the scalar develops a nonzero expectation value. (The field which 
we call $ later in this paper could be thought of as such a condensate.) 

4.2 Real u 

The observant reader will have noted that while both gravity arguments above rule out 
quasinormal modes with Imw > 0, they do not rule out quasinormal modes with Imu = 0. 
In Section |3[ however, we required the Green's function to be holomorphic also for Imcj = 0. 
In this section, in addition to ruling out quasinormal modes with real w / 0, we study some 
other properties of the Green's functions for real oj. 

Consider the quantity W = A' A — A' A introduced above, which for real oj can be as- 



sociated with the Wronskian. As discussed below (4.4), a necessary condition for A to be a 
quasinormal mode is that r n FW must vanish in the limit r — > oo. Recall that the Wronskian 
satisfies the differential equation 

W , (r)+T 1 (r)W(r) = 0, (4.6) 

which for T\ = (F' J 'F + n/r) is solved by W oc r~ n /F(r) up to an r-independent constant. 
Assuming A ~ (r — rh)~ %UJ l Fh near the horizon, we find 

r n F(A'A - A' A) = -iurl . (4.7) 

Thus the Wronskian never vanishes unless oj = (or in the limit r^ — > 0). We conclude there 
are no quasinormal modes on the real axis away from the origin for non-zero r^, and hence 
no poles in Gr(oj, k) for real 



9 We suspect that Y < also for a third example of instability: the real scalar fluctuation of the probe 
brane embedding in the D3/D7 system [Q. Based on numerical evidence, we know the Schrodinger potential 
of the scalar develops a well. 
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To understand the positivity properties of the spectral function, we use the Wronskian 
to show oj p(oj, k) > (for real oj). The idea is very simple. Up to an overall normalization, 
we can write the spectral density as 

Im G R (u)) ~ ^ Inn Fr n " t) ■ ( 48 ) 

Zl r->oo y A A J 

As A and A are linearly independent solutions, this expression is proportional to the Wron- 
skian W for the differential equation ( |2.7j ) . Thus the spectral density can be written in a way 
that makes its positivity properties manifest: 

TmG R (u)~% lim ^ , (4.9) 
2 r-^oo \A\^ 

using equation ( [4.71 ). 

To understand the reality property Gr(oj, k)* = Gr(—oj, —k) (for real oj), note first that 
( |2.7D is a real differential equation that depends on oj 2 and k 2 , not on oj or by themselves. 



The solution A is complex only because of the near horizon boundary condition (^2). Thus 
A, up to an overall constant phase factor, is a function of ioj with real coefficients. A is also 
an even function of k. Given the recipe for constructing Gr(oj, k) from A, the result follows. 

Note there exists a physical interpretation of the quantity fl4.7|) . 10 The action fl4.3|) 
contains a conserved Noether charge Q associated with phase rotations of A and A. The 
left hand side of ( |4.7| ) is equal to this Q. Thus ingoing boundary conditions at the horizon 
force Q ^ which in turn guarantees that there are no quasinormal modes for real oj 7^ 0. 
In the case where A is a scalar field, this charge Q is related to the radial component of a 
charge current J r . A slightly more elaborate argument relates (4.7) to current conservation 



V a J a = 0: J* is independent of t which in turn forces V r J r = -j^d r yj—g J r = 0. The left 



hand side of ( |4.7D is \J—gJ r up to normalization. The sign of Im Gr{uj) is related to the sign 
of the charge current J r . 

The absence of poles for real oj may seem surprising in view of some results in the liter- 
ature. For example, there are Dp/Dq brane systems in black hole backgrounds that support 



normal modes, i.e. poles at real oj at nonzero temperature [24, 25 1. Despite the nonzero tem- 
perature, these systems still have = in (I4.7I ). 11 Another distinct class of examples, where 
r/j = coincides with T = 0, has normal modes only in the limit T — > 0. At any T / 0, the 
modes have at least a small imaginary part. The current-current correlation function for the 
holographic superconductor with a charged scalar at the Breitenlohner-Freedman bound is in 
this class |29|]. Such modes were also observed in the low temperature limit of current-current 



correlation functions for the D3/D7 system, at nonzero charge densities 3C, 31, 32]. 



10 Ref. was the first to notice the importance of this quantity. 

n In the D3/D7 brane system with d = 4 at vanishing charge density and at sufficiently low temperature, 
there are (scalar and vector) fluctuations with modes having real frequencies. At low T 7^ the brane ends 
outside the black hole as opposed to falling into it. Thus the metric induced on the brane has Th — (Minkowski 
embeddings) [El E7I E| . 
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Even if we were to allow such a loophole of poles at real u 7^ here, given positivity of 
the spectral function uj p(u>, k) > 0, provided the poles are single poles, they will not spoil the 
sum rule. The single poles at oj ^ must have purely real residues. 12 Regularizing with an 
ie prescription, these poles in the real part of Gji(uj,k) introduce Dirac delta functions into 
the spectral density which we can integrate over. 

4.3 Contraction Map for Bosonic Correlators 

To study the large oj asymptotics of the Green's functions from gravity, we find it convenient 
first to convert Q2.7|) to an equivalent, nonlinear first order differential equation. If we have 
functions P(r),Q(r), then we can consider the quantity 

Q'A -QA' , 



(Usually we want P and Q to be approximations of solutions of (|2.7|).) Then 

,_ (sP-Q)[sV(P)-V(Q)} 

PQ'-P'Q • [ ' 



Also, we have the identity 

A' P's - Q' 



(4.12) 



A Ps-Q 
We define u = f r °° ^- , assuming that 

F ~ r 2 ( 1 - pr~ d + o{r- d ) ) (4.13) 



at large r and F ~ (r — r^) for r ~ r^. 13 Having set L = 2k 2 = 1, p here can be related to 
the on-shell value of the action and hence to the pressure of the field theory. For a review of 
this claim, see Appendix ||. We choose 



11 



1/2 



,1/2 



U 



Note that 



11 



W _W + _e\_r, (416) 



Q 2 \ 2r 2 2u 2 F 2 rF J F 



12 If the residue at oj = c has a nonzero imaginary part b, then p ~ b/{ui — c). The quantity u)p will change 
sign as uj passes through the pole. 



13 Equation (4.13) does not hold for all AdS theories. In general it will only hold if there are no fields with 
scaling dimension d/2 or lower. 

14 Equating the r.h.s. of equation (4.16) to zero determines those Y for which P and Q are exact solutions. 



This observation may be useful for constructing potentials with exact solutions. 
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and that the Wronskian has the form 

PQ> -P'Q = -JL^. (4.17) 

Thus the equation for s becomes 

s ' = £{s) = iS? { sH (n + m^ ~ H (n+i)/2( UJU )) 2 y (^) 

where we have defined 

(n (2-n n + 2 F' \ Y\ , lA . 

We will call the expression multiplied by n/2 above the pressure term because it depends on 
p through F. 

We will now derive an asymptotic expansion for s(r, u) in the limit uj — > oo, for Imw > 0. 
We introduce the mapping 

Z{s) = f £{s{r'))dr' . (4.20) 

We will assume that Imw > and that y(u) = 0(u 2A_1 ) as u — > 0. For the current-current 
correlator, we will assume A > (d — 2)/2 (the unitarity bound for scalar operators). For 
the stress-tensor correlator, we will assume that 2A > d. (This latter bound forces us to 
consider sum rules for the stress tensor with k = 0.) Given these assumptions, we will 
show in Appendix CT that this map ( [4.20| ) is a contraction mapping with contraction factor 



0(|o;| _1 ). Therefore, the error in s decreases by a factor of each time we iterate X. 

Generally, A is the scaling dimension of a field $ that is coupled to A. (Because of 
Lorentz and gauge symmetry constraints, both gauge and gravity fields typically couple to 
quadratic functions of fields, hence the factor of 2.) We are particularly interested in the 
effect of a nonzero pressure p, for which the effective A = d/2. (In the stress tensor case 
n = d — 1, the leading contribution from the pressure vanishes and A > d/2.) Although our 
emphasis on the gravity side of the duality obscures the point, one should be able to think 
of these couplings in terms of an operator product expansion (OPE) on the field theory side. 
The effect of p on the large frequency behavior of the correlation function comes from the 
presence of (T^ u ) in an OPE of the current-current or stress tensor-stress tensor correlation 
functions. Similarly, the effect of a quadratic coupling to a field <3? of dimension A corresponds 
to the presence of {OO) in the OPE (where $ is dual to O). 

4.4 Explicit Computation of Correlators 

Now, we actually compute the asymptotics of the correlators in question. 

The first observation is that the large to behavior of s is dominated by the large r region 
of the differential equation. This observation is intuitive from the AdS/CFT correspondence 
because large uj behavior corresponds to the UV of the field theory, and the UV of the field 
theory is precisely this large r region. From a mathematical point of view, we notice that 
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H^ n+l y 2 (oju) oscillates very rapidly when uj is large. If we expand a rapidly oscillating integral 
Jo°° f( u )9{ ljJU ) < ^ u as an asymptotic series in powers of uj, and / and g are C°° functions, then 
only their behavior near (the AdS boundary) and oo (the black hole horizon) is important. 
(One way to see this is to integrate by parts repeatedly.) In our case, / = plus some terms 
that are multiplied by F and therefore decrease exponentially for large u. So all derivatives 
of / go to zero at the horizon. 

In the proof of the contraction mapping above, we neglected suppression due to these 
oscillations. Typically, we expect that s will receive contributions only from the large r region 
of the integral and so be of order uj~ 2A . If the pressure term dominates, then 2 A = d while if 
the Y term dominates, then A will typically be the dimension of a scalar operator contributing 
to Y, for example the $ in flA.lOp . 

We will consider the current correlators in d = 3 and d = 4 and then the stress tensor 
correlator in general d. We begin with the current two point function in d = 3. From 
combining (2.11) and ( 4.12| ), we have that 



Gf(oj,0) = iu±-^ s , (4.21) 

(setting e = 1). Note that when d = 3, then n = and there is no contribution to s from the 
pressure term in ( [4.18 ). We have from the preceding discussion that a nonzero Y will induce 



an s ~ w _2A . We expect that the dimension should exceed the unitarity bound 2 A > 1. 
Thus, we expect the following variant of the Ferrell-Glover-Tinkham sum rule to hold for a 
generic choice of Y: 

o poo J, . 

Gf (0,0) = - / -(ImG?(,,0)- W ) . (4.22) 
n Jo uj 

For the current two point function in d = 4, we find that 

G|>,0) = uj 2 (ln(iA/ W ) - . (4.23) 

The contribution to s from the pressure term should scale as 1 /uj 4 while that from Y should 
scale as uj~ 2 ^ where now from unitarity A > 1. Thus a Ferrell-Glover-Tinkham sum rule of 
the form 

G£(0,0) = ~ jf ^ (**G%( U ,0) - (4.24) 

should hold for generic choice of Y . 

We can determine the leading contribution to s from the pressure term for this d = 4 case. 
If F = r 2 — pr 2 ~ d + ... for some constant p, then we have jp^z + -^i — = (— 4p/5)n 4 + .... 
Then we can integrate 

H^\uju) 2 — = — pcj- 4 . (4.25) 



o 



4 U 157T 



(The integral diverges when uj is real, but we may analytically continue from the upper half 
plane, where it does converge.) So the leading correction to G^(uj, 0) from the pressure term 
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is (— i-iru) 2 ) (jff^pw 4 ) = (8p/15)a; 2 . This term agrees with an analytically solvable case 
discussed in Appendix A.l| . See in particular ( A.33| ). By dimensional analysis, this correction 



is of the form one would expect if the first nontrivial term in the OPE of J x with itself 
contains the stress-tensor. 

For the stress tensor correlator, we are interested in the component G x R v,xy {uj, 0) = Gr(uj). 
Combining the results of Appendix [B] and ( |4.12| ), we find that 

G*M = -p-f W-^J, (4.26) 

G R (u>) = -p - f + \^ ln(<A/ W ) - , (4.27) 

in d = 3 and d = 4 respectively. Similarly to the current-current correlation function, the 
presence of p in the expressions above is consistent with what one would expect if the OPE 
of the stress-tensor with itself contains the stress-tensor. To eliminate the power law growth 
with respect to co and construct a sum rule, we consider the regulated Green's function 

5G R (u>) = G R (u, T) - G R (u, 0) , (4.28) 

where p and s will generically depend on temperature T. The leading p in Gr(uj) comes from 
a contact term that we determined in Appendix [B[ The second, p dependent term comes 
from a correction to Q. This correction can be generalized to arbitrary d. We have to look 
at the large r limit of Q: 

Q = (w)' 1 -^ 2 / 2 ^^™) . (4.29) 
and in particular we have (r^M/ 2 = 1 - j£j^r- d + 0( r - d - 2 ). So Q - Q *=° - z^ll^-d 



2{d+i) 1 ^^V' )■ ^ Qp = - 2 (d+l) 



r 



and 6G R (u) =-p- p{d - \)d/2(d + 1) + 0{s) 

To understand the correction from s, note that having set k = (to be consistent with 
the bounds on A and have a contraction map), we find Y = in ( 4.18[ ). As mentioned above, 



in the special case n = d — 1 the leading contribution to s from the pressure term vanishes. 
We would need to specify more terms in a power series expansion for F to compute s reliably, 
but we know that s ~ u>~ 2A where A > d/2. Thus s will lead to corrections of Gr(uj) that 
vanish in the limit cj — > oo. We find 

l^5G R (u) = -p-^-^. (4.30) 
^->oo 2{d + 1) 

To understand the sum rule, we also need the behavior of G R y ' xy (0, 0). When w = 0, the 



equation (2/7) is easy to solve, A = const. Thus from Appendix ||, G x R v,xy (0, 0) = —p, and we 
find the sum rule 

- I - W = { 4^P ■ (4-31) 
vr J u> FV ; 2(d + 1) v ' 

In d = 4 we reproduce the result of Romatschke and Son Q: ^ J ^-5p(uj) = |p = |e. (In the 
last equality, we used the fact that for a conformal field theory, the stress-tensor is traceless 



and so the energy density e = 3p.) Note we expect this result ( 4.31 ) to hold in arbitrary d 
for a generic F satisfying our assumptions. 
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5. Fermions 



We now repeat our analysis carried out in the previous section for the case of Green's functions 
of fermionic operators. Again, it has to be shown that the quasinormal modes of (now 
fermionic) bulk fields exclusively lie in the lower half of the complex plane and that the 
spectral function falls off fast enough at large frequencies. Our strategy is to decompose the 
Dirac equation into a set of identical differential equations involving two-component spinors. 
This decomposition allows us to study sum rules in arbitrary dimensions in a simple way. 

5.1 Fermionic Green's functions in arbitrary dimensions 

Consider a fermionic operator 0<i> in a field theory dual to a spinor ^ in gravity with the 
action 

5* = -i / d d+1 x ^j^(~f a D a - , (5.1) 



where \t is a Dirac spinor of charge q, \P = ^'7*. $ is a real scalar field which allows 
us to study fermionic fluctuations in presence of a scalar condensate, e.g. in a holographic 
superconductor. 15 We do not separate <3? from a possible fermion mass m. When we eventually 
take the metric to be of the form ( |2.5| ), we will assume that $ goes to zero at the boundary 
as rah + <I>o r ~ A for some A > 1 (which is in general different from the bound on A in the 
previous section) and that 3> is finite at the black hole horizon. For the sake of generality we 
keep the fermion mass m 7^ for now, but will set it to zero soon below. Furthermore, we 
assume that At is the only nonzero component of A. The covariant derivative is 



D a = d a + l -u a fi_cl- ~ iqA a , (5.2) 

where oj is the spin connection and A a is a gauge field. For a diagonal metric that depends 
only on the radial coordinate r, the spin connection can be written 



r]^\fF ¥ d r \J\g^u\ • (5.3) 



The gamma matrices 7° satisfy the Clifford algebra {7", j b } = 2g ab . Vielbein indices are 
underlined, and the generators of the Lorentz group are a— = ^[7-, 7-] = ij—. The dimension 



2 

of is d/2 ± mL, where A^ > (d — l)/2 satisfies the unitarity bound [34] 



Dirac equation in arbitrary dimensions Assuming the metric is diagonal and depends 
only on the radial coordinate r, it is convenient to rewrite the equations of motion in terms 
of a rescaled (2 n )-component spinor tp = (—g ■ g Tr ) 1 ^'^ where n = [(d + l)/2j. Given the 
above assumptions on the form of the metric, the equation of motion for if) is 

[-/ a (id a + A a )-i<S>]iP = , (5.4) 

setting q = 1. 



3 This kind of coupling along with others, and their various implications, were discussed in 
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We suppose that the fermion has some momentum k. Because of rotational symmetry, 
we may assume without loss of generality that the momentum is in the x direction. The 
equation of motion for ip is then 



sf^frf-iu) + A t ) + i^l-dr 



l k -i<& 



ij) = 0. 



(5.5) 



We argue that ( |5.5| ) can be decomposed into a set of identical differential equations 
acting on two-component spinors. Consider a Euclidean Clifford algebra, {j a , ; y b } = S ab in D 
dimensions. In the case D = 2, we take 7 = 03 and 7 1 = 02 the Pauli spin matrices. Given a 
Clifford algebra j a in D = 2n dimensions, a Clifford algebra r a in D + 2 = 2n + 2 dimensions 
is 



id< 



a 3 , i~ = id®ff 2 , r a = 7 a ~ 2 ® ax . (5.6) 

For odd dimensions, in the usual way we identify ^ 2n+1 with the product of the other gamma 
matrices (up to a factor of i). Given this construction, we can think of 7 and 7 1 as block 
diagonal, consisting of 03 or 02 matrices along the diagonal. Additionally, 7 2 (for D > 2) 
is block diagonal, consisting of o\ matrices along the diagonal of alternating sign. The 
Lorentzian Clifford algebra can be recovered by multiplying one of the 7° by i. We choose 
7- = ry 1 , 7- = 7 and 7- = 7 2 . With these choices, ( |5.5| ) in an arbitrary number of dimensions 
d reduces to two decoupled equations 



y / -g tt ia 2 (uj + A t ) + i\fF F CF3d r + {-l) a ^/g^aik - i $ 



(5.7) 



where i\) a with a = 1,2,... ,2 ra_1 are the 2( n_1 ) two-component spinors appearing in the 
2( n_1 ) spinor equations, respectively, with momentum alternating between +k and —k. There- 
fore the second block of equations, i.e. the one for ip2, is related to the first block by k — > —k. 
A third block (if present) is identical to the first, the fourth block to the second, and so on. 

Fermion operator Green's functions We are interested in the retarded Green's function 
for the boundary fermion operator O^. The prescription we use to study this fermion corre- 

. In 



lation function closely follows that of [35, |36], 37, f38[ which is based on the work of j 
what follows we will summarize this prescription. Consider a solution to the Dirac equation 
with infalling boundary conditions at the black hole horizon. Before setting ml = below, 
let us consider small masses mL < 1/2 now. Decomposing ip into eigenvectors of 7 r , for 
mL < 1/2 we find the AdS-boundary behavior 




a a r 
b n r 



•-mL _|_ Q(ymL—\ 



-mL 



+ Oir 



-mL—1 mL—1 
1 ' 



a 



1, 2, 2 



(n-l) 



(5.8) 



We are going to work in the massless limit from now on, so mL 
components scale with the same power of r near the boundary. 16 



and the two spinor 



Note that the boundary expansion ( p.fij) can involve logarithmic sources for half integer masses. Further- 
more, for mL > 1/2 the term involving b a r~ mL would become subleading compared to 0(r mL ~ 1 ). 
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Note that for mL = two quantizations are possible, choosing either a a or b a as the 
sources. Here we choose to identify the number a a with the source for the operator O a q,, 
while b a = (O a ^). Since the Dirac equation is linear, a a (uj,k) will be linearly related to 
bp(cj,k). In order to make this relation explicit, we define the 2 n ~ 1 component spinors Tp±. 
At the boundary, ip + asymptotes to the sources, ip- to the vevs: 



lim ipj r = a ; lim ip- = b . 

r— J-oo r— >oo 

The linear relation between vevs and sources can be written 



(5.9) 



b(uj, k) = S a(oj, k) . 



(5.10) 



The implicit prescription to compute the retarded Green's function is then [35]: 



-iG R -fl dy Ma 



(5.11) 



Note that 7bd y * s the field theory gamma matrix. M is a change of basis matrix that allows for 
a and b to transform in different representations of the Clifford algebra. In odd dimensional 
AdS spaces, 7^ dy = i.e. 7- restricted to the positive eigenspace of 7-, and M = id. In 
this case Gr is defined only on the negative eigenspace of 7- and 7- itself is reinterpreted 
as the boundary gamma matrix which determines the chirality of the field theory spinors. 
With our boundary conditions, the expectation values b are negative chirality Weyl spinors. 
Choosing the other quantization would allow us to study the opposite chirality. 

For even dimensional AdS spaces where d+1 = 2n, we need to construct the field theory 
gamma matrices. One simple choice is to let b a transform under 



Tbdy 
Tbdy 



7 C 7bdy 



such that a a is related to b a via M 
and 7— are compatible. 



7h>dy 



We rewrite the correlator recipe (5.11) 

b, 



With this choice, the Lorentz generators 7^ 



^, a €{1,2, ...,2(- 1 )}. 



(5.12) 



Note that in the present setup the matrix-valued Green's function is diagonal, because the 
spinors tp a decouple due to the block diagonal structure of the Dirac equation discussed above. 
We now analyse the pole structure of the fermionic Green's functions from the gravity point 
of view in analogy to the bosonic case in the previous section. 



5.2 Imw > 



Analogous to the bosonic case discussed in Section 4.1 we can again argue for the absence of 



branch cuts as shown in Appendix D.2. Thus the only non-analyticities to consider are again 
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the poles of Gr which now correspond to quasinormal modes of tp. As was the case with the 
bosons, there are no poles in the upper half of the complex plane (Imw > 0). Quasinormal 
modes are solutions to ( |5.5| ) that are normalizable at the boundary (a a = 0) and ingoing at 
the horizon. In the upper half plane, the solutions are also normalizable at the horizon. So 
quasinormal modes in the upper half plane are eigenfunctions of the Dirac operator of (5.5) 



with eigenvalue uj. But in the upper half plane the Dirac operator is self-adjoint and can have 
only real eigenvalues. 

As mentioned in Section because is not Hermitian, the retarded Green's function 
is not expected to have a reality property of the form G°^ {uj,k)* = G°^ {— uj, — k). For our 



fermion, however, the Dirac equation (|5.7l ) has a symmetry under complex conjugation plus 
uj — > —uj, k — > —k and At — > —At. Thus the diagonal components of the Green's functions 
are expected to have the symmetry Gr(uj, k, p)* = Gr(—uj, —k, —fi) where p is the boundary 
value of A t . 

5.3 Real uj 

As was true for the bosonic Green's functions, we now show that there are no poles on the 
real axis, i. e. no (quasi)normal modes with real frequencies. The fermion spectral function 
can be written in terms of bulk field components 



p(co,k) a a = ^[G Ra a -(G Ra ^ 



lim — 

r— >oo 2 



oj,k£K. (5.13) 



The spinor components ip a ± are chosen such that near the horizon they behave as ip a ± ~ 
( r — r h)~ lul l Fh [ca,± + 0(r — rfr)]. p has poles if ip a ,+ and ip^ + vanish simultaneously. The 
quantity Wf = (ipa-ip^ + + ip a)+ ip^ _) is the fermionic analog of the Wronskian considered 
in Section |], and has to vanish at a pole. It follows from the equations of motion that Wf 
is independent of the radial coordinate, d r Wf = (see e.g. [39, ^]). 17 Exactly at the 
horizon we find for the diagonal entries 

lim W f = -2\c a+ \ 2 = -2|c Q _| 2 , (5.14) 

because, at uj ^ 0, the leading coefficients are fixed by the infalling boundary condition to 
satisfy c Qj _ = — c Qj +. Thus there are never poles in the fermionic Green's function along the 
real axis away from uj = (c Qi ± = is the trivial case ip a ,±( r ) = 0). In contrast, at uj = the 
infalling boundary condition and equations of motion no longer imply a relation between the 
leading coefficients c a ±. In this case, without loss of generality one of these coefficients can 
be chosen to be real, the other imaginary, and the sum Wf = [c* _c a) + + c a -C* a + ] vanishes. 
For example, in the zero temperature limit a fermionic (quasi)particle pole sits at uj = at 
the Fermi momentum k = kF 



17 Similar to the bosonic case, this quantity can be related to the radial component of a conserved current 
jr _ qyyijr = Wflsf—g- That d r Wf = follows from the fact that V a J a = and the observation that the 
components of J M depend only on r. 
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Positivity of the spectral function can be made manifest using equation ( fug) Q. Con- 
sider that p a a = — lim^oo Wf /(2\ip aj+ \ 2 ), which implies for the diagonal elements that 

Pa a = lim |Cq '+ I ' >0. (5.15) 
r-^oo \ip a:+ (r)\ 2 



Note, that the fermion spectral function (5.15) does not contain any explicit factor of oj, in 
contrast to the bosonic one (|4.9|). Therefore, it is non-negative also for negative frequencies. 18 



5.4 Asymptotic Series for the Correlator 

Again, we want to determine the large frequency behavior of the Green's function for Im ui > 0. 
In order for a sum rule to be valid, we need Green's functions that fall off fast enough at large 
frequencies. In order to examine this large frequency behavior, we are going to introduce a 
map Let us first define two helpful abbreviations 



and 



dr'\/-g u g r r , (5.16) 

poo 

/ dr'y/-gttg rr A t . (5.17) 

J r 

While the method that we used for bosons can also be applied to fermions, 19 we choose 
to use a different method that can be generalized more readily to spinors with more than two 
components. We define the mapping 

with 



00 / / r,XX 



Jty) = ie^H^u+viu)) j dule -iy^t {ulu > +v) ^r ^9^_ kry x + _J^_j ^ ( 5 1Q ) 

where satisfies 

7 ^(°) = _^(°), d u ft 0) = 0. (5.20) 



Then a fixed point of X is a solution to ( |5.5| ) with infalling boundary conditions. In other 
words = ip is equivalent to the equation of motion (5.5), as can be shown by taking the 
derivative with respect to u on both sides of = ip. 



18 Similar to the bosonic case, the sign of the spectral density is related to the direction of the radial current 
J r mentioned in footnote [nj. 

19 Suppose A, P, Q are functions that take values in a two-dimensional vector space, and A satisfies T>(A) = 
A' + TA= 0. We define s = Then s' = ^ [(P A V{P))s 2 — (PA V{Q) + Q A V{P))s + (Q A V(Q))} . 

a\ ( P 

The second order system A" + T\A' + To A = is just a special case with A — > [ Al ) , P 
Q ,YandX ^° ^ 



A' \ P' 



T T 
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We can write the first term of Q5.18 ) as 

^(°) = 1(0) = e-iM+^o) . (5.21) 
This is a solution to the equation of motion (|5.5|) for k = and <3> = 0. Therefore it can also 



be considered an approximate solution to the full equation (5.5). We will show in Appendix 
|C.2| that the map I(ip) can be iterated, using 

^(0) 

as an initial approximate solution, to find 

the exact solution. 

Now we actually compute the first correction to the fermion correlator in the large uj 
limit. The leading term comes from a ratio constructed from the components of ip^ and 
ipt 1 ) = I^ijj^) — i/>( ). We compute V^ 1 ^ and obtain 

ibM = ie i(^+v) f°° du , e -2i(uu'+v) r ( [g^_ kl x + _J^ = \ ^(0)_ ,5 22) 

Ju \\ -9 tt V^g* J 

Once again we have an integral that oscillates rapidly for large uj. As long as all derivatives 
of A M 7~ (\/g xx k"f- + i$>) go to zero at the horizon, the asymptotic series depends only on 

the behavior near the boundary. We have \/—g tt ~ \/g xx ~ u. If we let At ~ /x, then v ~ /xu, 
so v — > at the boundary. We also let <3> ~ <3?o^ A with A > 1. 
The leading correction from k is 

(5-23) 

The leading correction from <3? is 

4\ =0 ~i-^r(A) 7 ^(°). (5.24) 

In our gamma matrix basis, we choose ^(°) = (1, 1, . . . , 1) <g> (-1, 1). We find then that 
^W| u=0 ~ (fc/2w)(l, -1,1,-1,... ,1,-1) ® (1,1) +o(l/w). Therefore at large u 

G R (u, k) a a =U±*+ o(l/u)j . (5.25) 

Note that if we consider instead the trace (provided d > 2), 

tiG R (u f k) = 2 n - 1 i + o(l/uj) . (5.26) 
Thus, we can write down a sum rule of the form 

/j 
— (trp(u,k) -2 n ~ l ) . (5.27) 
7T 

The left hand side will typically vanish provided tr Gr(uj, k) is finite at the origin. This sum 
rule is similar in spirit to but different in detail from a fermionic sum rule that arises in the 
context of ARPES (Angle Resolved Photo Emission Spectroscopy) [^Tj. This ARPES sum 
rule is j dcjA(u},\c) = 1, with the spectral function A(oj,"k). 
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6. Discussion 



We have seen how some well established properties of Green's function in field theory are 
realized in a dual gravitational description through the AdS/CFT correspondence. While 
our results are not completely general, they do apply to a class of correlation functions 
involving scalars, conserved currents, the stress tensor, and massless charged fermions. Some 
of these properties were easy to establish and the results were partially known in the literature 
before. For example the positivity properties of the spectral function, ujp(uj, k) > for bosons 
and p(u), k) > for fermions, follow from the horizon boundary conditions and a conserved 
charge. That Gn(u>,k)* = w, — /c) for Hermitian operators turns out to be an almost 

trivial feature of the differential equation and horizon boundary condition governing the 
Green's function (and similarly Gn(Lu,k, /z)* = Gn{—uj,—k,—p) for fermions). Other results 
were more significant. For the bosons and fermions, we established necessary conditions for 
Gr(lo, k) to be holomorphic in the upper half of the complex frequency plane including real 
cu 7^ 0. We developed a contraction map for the bosons and a contraction-like map for the 
fermions that allowed us to determine the large frequency behavior of Gr(oj, 0). We suspect 
that these contraction maps may be useful in other contexts to study the large frequency 
behavior of Green's functions from gravity. 



We were able to rederive some known current and stress tensor sum rules ||, 11, 12] but 
from a different viewpoint and in a more general context. From the results of this paper, we 
see that these sum rules hold not just for black D3-brane backgrounds || |l^| and not just 
for the holographic superconductor in AdSi with a scalar operator of conformal dimension 
one or two [11], but more generally, in other space-time dimensions and in cases where the 
backgrounds may be deformed by expectation values of arbitrary operators (though typically 
some lower bound was placed on the conformal dimension of these operators). 

One interesting observation that follows from generalizing the sum rule of || to arbitrary 
dimension is that fl4.3l| ) for the stress-tensor correlator G xy ' xy (oo, 0) can be related to a similar 
sum rule for Chamblin-Reall backgrounds |J. The sum rule of ||] can be written in our 
notation as 

fir*? («> 

As discussed by |42j[ , these five dimensional Chamblin-Reall backgrounds can be realized as 
dimensional reductions of asymptotically AdSd+i space-times on a d — 4 dimensional torus. 
The speed of sound v s and temperature T of the Chamblin-Reall background are identified 
with those of the AdSd+i space-time, while the entropy density s and pressure p are rescaled 
by the volume of the torus. For conformal field theories dual to these AdSd+i space-times, 
we expect sT = pd and = l/(d— 1). Making these substitutions converts (O) into ( |4.31 ), 
As pointed out by ||43|| , this trick of dimensional reduction also works for constructing 
black Dp brane solutions from higher dimensional AdSd+i space-times. We find that black Dl 
and D4 brane solutions come from reducing AdS^ and AdSj space-times on a circle, as one 
might expect given the relation between M-theory and type IIA string theory. More formally, 
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DO branes and D2 branes come from AdSd+i where d = 14/5 and d = 10/3 respectively. 
Thus we find predictions of Romatschke-Son type sum rules for black DO, Dl, D2, and D4 
branes. In general, Dp-branes are related to dimensional reductions of AdSd+i space-times 
withd = 2(p-7)/(p-5) @. 

In the appendices, we wrote down some mildly interesting sum rules that to our knowledge 
are new. We noted that the current correlation function in d = 2 has a pole at oj = that 
makes the integral in the naive sum rule ( |3,7j ) divergent. One possible construction that evades 
this problem is flA.26 ). Another more standard sum rule ( A.27| ) involves regulating the pole 
by taking derivatives. We also wrote down some new sum rules ( A.37| ) for the current-current 
correlation function in the d = 4 case, i.e. the case of N = 4 super Yang-Mills at large N and 
large 't Hooft coupling. We hope our techniques can be used and generalized to find more 
new sum rules from gravity. 

We did not study the full current-current or stress tensor-stress tensor correlation func- 
tions, but only a single component of these multi-index objects. Our fermionic Green's 
function turned out to be diagonal in spinor indices. One would of course like to investigate 
the full current-current and stress tensor-stress tensor correlation functions, and also consider 
more generic cases where the fermionic Green's function is not diagonal. Such cases would 
involve solving more general coupled systems of ordinary differential equations in gravity. 
Two distinct systematic approaches to such coupled systems were developed in 46]. 

We leave the reader with the most interesting unanswered question of the paper: What 
is the physical meaning of the sum rule (3/7) in gravity? Is there a dual gravitational inter- 
pretation of the sum rule itself? 
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A. The Holographic Conductivity Equation 



In this appendix, we want to demonstrate how an equation of the form ( \z. (\) arises m com- 
puting a current-current correlation function from AdS/CFT. We also review some simple 
cases where closed form expressions are available for these correlation functions. 
Consider the action 



S 



d+l, 



-9 



1 

2^2 



R + 



d(d - 1) 
I? 



1 

4? 



F ab F' 



ab 



\Da*\ 



m 2 |<3?| 2 
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where D a = V a — iA a . We look for a solution to the equations of motion with the ansatz 
A = A t (r) dt + A x (t, r) dx, $ = 3>(r) (with $ a real function) and 



ds 2 = -f(r)e-*Vdt 2 + -jo + ttt + 2 fe (i, r) dt . 

L 2 /(r) 



(A.l) 



We will treat the components A x and gt x as perturbations and only consider their linearized 
equations of motion. As given in ref. |i"4j| , we find that the background satisfies the equations 
of motion 



d>'2 + 'LllJL + J_^L + (d-l)(d-2) _ d(d-l) | 



/ 2e 2 / 



2k 2 t 2 
+ 

d-1 x' 



e*A 2 § 



2n 2 L 2 f 

+ 



2 2k 2 



_ 2 re*A 2 $ 2 



4' + 



/ 2 

2e 2 A 4 $ 2 

7 









$>" + 



h ^ - ^- $ + 



/ 



; *A 2 



f 



m . 

-=- $ = 



We impose the boundary conditions 



(A.2) 
(A.3) 
(A.4) 
(A.5) 



,-A 



' X f = h 

L 2 



$ = Q A r A + . . . , A t = /it + . . . 



(A.6) 
(A.7) 



as r — > oo where A > (d — 2)/2 satisfies the unitarity bound. We also want to assume the 
existence of a horizon at a radius r = where /(r^) = and A^r/J = vanish linearly 
with r. The other functions x and $ are finite at the horizon. There is a relation between 
A' t {r h ) = A h , f(r h ) = f h , x(r h ) = Xh and *(r fc ) = 



'7, 



L 2 d- 1 



A 2 ^" 



+ 2m 2 $| 



(A. 



Given such a black hole background, the linear fluctuations satisfy the differential equa- 



tions 



9tx 9tx H y t * 



^_2^Af 
n e 2 / 1 ' 



4,+ 



2e 2 $ 2 



,/ 2 e 2 / J f 
This equation ( A.lOj ) is precisely of the general form fl2.7| ). 



A, 



, 



. 



(A.9) 



(A.10) 
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In the following subsections, we review the cases d = 2, 3 and 4 when = At = x = 0, and 
/ = r 2 (l — (r^/r) d )/L 2 , and the equation ( |A.1Q ) is exactly solvable. One can locate the poles 
of the correlation function, deduce the large u behavior, and verify the sum rules analytically. 
The sum rules for the corresponding correlation function G^{uj) are often rephrased in terms 
of the charge conductivity a(u)) = G x ^{oj)/iu. 

A.l Current-Current Correlators in Thermal Backgrounds 

We review some known and partially known results for current-current correlators calculated 
from dual AdS-Schwarzschild backgrounds. These are correlation functions in a thermal bath 
that preserves rotational symmetry. The first part of our treatment follows [48|. 

Gauge invariance puts some constraints on the form of current-current correlation func- 
tions. We expect that 

G^(u, k) = P^nV, k 2 ) + P^U T (u, k 2 ) , (A.ll) 

where 

i£ = 0, P£ = 0, P? 3 =5 t3 - k -§, (A.12) 

k z 

huh,, 



P^ v — Pfii/ P^ u , P/iu — f][i,v ^ 2 - (A. 13) 

These correlation functions satisfy the gauge symmetry Ward identities because k^P^ v = = 

Our d + 1 dimensional black brane metric is 

ds 2 = ^ {-f{z)dt 2 + dx 2 + j^) , (A.14) 

where f(z) = 1 — (z/zh) d - Consider the action for a Maxwell field in this background: 

Sem = -\J d d+1 x^F ab F ab . (A. 15) 

We would like to consider fluctuations with a e ~ lujt + lkx dependence. There are two types, 
longitudinal and transverse. By a gauge choice, we set A z = 0. The transverse fluctuations 
involve A a where a ^ t, x, or z, while the longitudinal fluctuations involve both At and A x . 
The transverse fluctuations satisfy the differential equation 

7' 3- cA u 2 - k 2 f 
J + ~) a + _ f 
The longitudinal fluctuations satisfy the three coupled differential equations 



<+( i 7 + —)K + ^^ L A a = 0. (A.16) 



ujA't + kfA' x = , (A.17) 

/' 3-d \ , lo 2 uk 
J + -) A x + T2 A x + j2 

K + 3 -^A' t -^At-^A x = 0. (A.19) 
z J J 



A " + I - + — ) A X + 72^- + 72 A t = » ( A - 18 ) 
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The third equation follows from the first two. There is an equivalent "gauge invariant" 
formulation of this system of equations constructed from E = kA t + ujA x : 



E" + 



0J 2 f 



+ 



3-d 



f(uj 2 — k 2 f) 



E' + 



uj- 



k 2 f 



P 



E = 



(A.20) 



Prom solutions to the differential equations ( [A, 16 ) and (A.20), we can construct n T and 
Tl L respectively. For a retarded Green's function, we choose ingoing boundary conditions at 
the event horizon: E ~ A a ~ (1 — z )-' lu>z h/d_ sma Ji z , A a and E have the expansion 

o(l + ...) + 6^- 2 (l + ...) , (A.21) 

and the corresponding retarded Green's function, up to contact terms, should be b/a for 
d > 2. For d = 2, the expansion will be a little different, b — a In z + . . . (see for example |]49| ). 

In the following subsections, we will work in units where z^ = 1. Note that the tempera- 
ture of the field theory is given by T = d/Airzfr. To restore units, we shift uj —> uiz^ = uid/AirT. 

d=2 

In d = 2 space-time dimensions, there are no transverse fluctuations, and we need solve only 
(A.20). The solution with the correct ingoing boundary conditions at the event horizon can 
be written in terms of hypergeometric functions 



E 



:i-z 



'l-z 2 )z 2 



or - k 2 



+ ioj — 1 x 



(A.22) 



x 2 Fi 



i(uj + k) iioj — k) o 
2 - v Z , 2 ~ ~ ~, 2-iuj;l-z 2 



+ 



2 — iui 



1 (iw - 1) 2 F 1 



(The other solution is the complex conjug ate.) 20 The near boundary expansion of this solution 
has the form E = b — alnz. The corresponding Green's function is then 



U L (oj,k) 



■iui 



2 - k 2 2^1 2 



i(k + uj) 



7 + In A , 



(A.23) 



where In A is a renormalization dependent contact term and ip(x) = T'(x)/T(x). Note that Tl L 

has poles in the lower half of the complex plane at positions uj = ±k—2ni where n = 0, 1, 2, 

The spectral function is 



p{uj,k) =Imil L = n - ( coth [—(uj + k)) + coth (—(uj - k) 



(A.24) 



Since uj p(uj , k) / (uj 2 — k 2 ) > 0, it follows that ujp(uj,k)P L is a positive definite matrix, as it 
should be. 



20 We would like to thank Jie Ren for discussions about this case. 
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There is a sum rule here of the form 21 

00 did ( , , , ,, . b — a\ it b , . 

— p(u)/a, - p(u)/b, + U-b-. A.26 

to \ oj J 2 a 

Thinking of a and b as shifting the temperature, we are regulating the sum rule by comparing 
the spectral density at two different temperatures and additionally subtracting off a divergent 
contribution at to = 0. Note that p(oo,0) = n/2. There are also a series of generalized sum 
rules of the form (p 

8 J 7TUT \ 

^ (4) (i) r duj ( , . 1 ^ - 1 - :! 

' 'p(^O)--- — + 



384 7 vra; 5 V w 12 720 

t// 6 )(1) f doj / . 1 ^ 2 w tt 4 cj 3 vr 6 w 5 
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7r 2 c<; 


Co' 


12 


1 


7r 2 c<; 


Co' 


12 


1 


7T 2 CU 


w 


12 



46,080 J ttlo 7 V ' w 12 720 30,240 



^ (2n) (l), /" / , n . 1 l^(™) 2fc B 



■1 



/" / , , 1 If iKujy K B 2k \ . . 



(4n)!! v ' 7 7rcj 2 ™+ 1 \ rv ' ' cj w ^ (2fc)! 

The subtractions here regulate the IR divergence at to = 0. 
d=3 



For (i > 2, an exact solution is available for ( A. 16| ) and ( |A.20| ) only in the case k = 0. Note 
that when k = 0, the two equations are identical and n L (u,0) = n T (cu,0). 
For d = 3, the authors of found the solution 

A a (z) = £(z) = exp Uj £ j^rj . (A.28) 

The corresponding Green's function is 

II(w, 0) = ip = ioj . (A.29) 

Note that the Green's function in this case has no poles. The spectral function satisfies the 
trivial sum rule 

- l) duj = . (A.30) 



21 This integral is an example of Frullani's integral (which in turn is an application of Fubini's theorem). 
Assuming that df/dx = g(x), then 

dx f(bx) - }{ax) = f 6 ^ f 00 ^ = _ _ (A.25) 







In our case, f(x) ~ coth(x) — 
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d=4 



For d = 4, the authors of [52] found the solution 

1 + 



1 



(A.31) 



2^1 



1 + « 1 + i !W 1 . 
1 OJ, OJ, 1 , - 1 

4 4 2 2 



(The other solution is again the complex conjugate.) The Green's function is 



n(w,o) 



-10J 



27 + In A' + V 



1 + t 



(A.32) 



4 J ' ' V 4 

This Green's function has poles in the lower half of the complex plane at positions oj = 

2n(±l — i) where n = 0, 1, 2, For comparison with a result in the body of the paper, we 

give the large oj expansion of H(oj, 0): 



n(w,0) = u} 2 ln(iA/u) + 



15w 2 



+ O(oj- 



The spectral density is 



ttoj 2 sinh(-7TLij/2) 



P 



An indefinite integral is known 
,2 



/ 



2(cosh(-7ru;/2) — cos(7rcj/2)) 



oj 



(A.33) 
(A.34) 



5^ = ^ + |(l-z)ln(l- e —^)/ 2 

+l L i 2 ( e -™(m/A _ l L i 2 (V™^ 



• =■(! • /) 111 (1 , ; 2 



■0/2 



From this integral, it is straightforward to verify the sum rule 



OO -y 



TTOJ 



duo = . 



(A.35) 



(A.36) 



r oj v 2 
This integral was studied numerically by |i"2|] . 

Like in the d = 2 case, there are some more interesting sum rules of the form ( |3.8| ) that 
one can write down: we have 

C(5) f duo ( ^ _ tt 2 oj 3 tt 4 oj 5 \ 
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uj + ojJ2(-^ 
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n 2 0J 2 (4k-2)\ (4fc)! 
The subtractions regulate the IR divergence at oj = 0. 



(A.37) 
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B. The Holographic T xy Two Point Function 



In this appendix, we review how an equation of the form ( \2.1\ ) arises in computing a stress- 
tensor two-point function from AdS/CFT. Consider the action 

S=±J ^xV=g (R + + ^ / Af^ (A" - ^1) , (B.l) 

where K = V^n^ (sum on only the gauge theory indices) is the trace of the extrinsic curvature 
and is a unit vector normal to the boundary and pointing toward larger r. We take the 
ansatz for the metric 

™2 J™2 „2 

ds 2 = -f(r)e-^dt 2 + -jdx 2 + i— + 2-^0(r, t, z) dy . (B.2) 

L z j(r) L z 

The tensor gf^ry is the induced metric on a constant r slice. 

We define two quantities, the bulk Lagrangian and the zeroth order Einstein tensor: 



(B.3) 

0=0 



The first observation is that 

£ _ (2 _ ^ )G ^ = _ i [^(a a 0)(9» - 49 a (^^»] 



4k 2 

a=t,r,2 

k z \ r y 

Note that this equality is only valid to order (f) 3 . From this observation, because of the (d(f)) 2 
term in the Lagrangian, we conclude that ( |2.7| ) does indeed govern the stress-tensor two-point 
function. Also, we find the onshell action reduces to 

S os = ± J d d x^g^ (K-*=±- n *(±- \<t>d r <^ ) . (B.5) 

With the addition of matter fields to the action, one should include the bulk stress tensor 
in making the subtraction ( |B.4 ). However, provided these matter fields have high enough 
scaling dimension, they will not affect ( |B.5| ) in the large r limit. 

Consider the simple case of a constant metric perturbation where eft = 4>o is a constant. 
We will assume that for large r, x ~ r~ 2A where A > d/2 and that / ~ (r 2 /L 2 )(l — 
2K 2 L d+1 p/r d + ...). In this case, the on-shell action evaluates to 

S os = j d d xp (l - ItfA . (B.6) 

We conclude two things. The constant p can be interpreted as the pressure and G x ^' xy (0, 0) = 
-p. 

To see what happens when cp is space-time dependent, let's consider the specific cases 
d = 3 and d = 4. Again, this analysis should generalize to the case where there are additional 
matter fields in the action, provided the matter fields have conformal scaling dimension A > 
d/2. 
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d=4 



Near the boundary, (j) has the expansion 



This expansion yields the on-shell action 

\u 2 - k 2 



16r 4 
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(B.7) 
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To regulate the ln(r) and r 2 divergences, we can add boundary counter-terms 

L 2 In A' 



Sct = i 2 / d4x V~^y ( n - 



— TV xv 'R, iXV j . 



where TZ^ and 7?. are the boundary Ricci tensor and Ricci scalar respectively. At 
order in 0, 



The end result is that 

<Sos + <Sct 

d=3 

Near the boundary, <f> has the expansion 
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(B.10) 
(B.ll) 
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This expansion yields the on-shell action 



/ 
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To regulate the r divergence, we can add a boundary counter-term 



Sct = Jk 2 j d4x V~9hryK , 



where 1Z is the boundary Ricci scalar. At leading order in 4>, 
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The end result is that 
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(B.14) 
(B.15) 
(B.16) 
(B.17) 
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C. Contraction Maps 

C.l Bosons 

We prove 

X{s)= r £(s(r'))dr' (C.l) 

Jr H 

is a contraction mapping. First, we restrict the set of s that X acts on to a domain D = {s : 
||s||oo < K} for some constant K where ||s||oo = su Pre[r H ,oo[ \e~ 2%u)U s{r)\ We will show first 
that I(s) G D. Then, we will show that for sufficiently large cj, X is a contraction mapping 
with contraction factor 0(|w| _1 ). More specifically, we show that 

|J(fl)-J(fl )| < A||*-ao||oo Vs,s GD, (C.2) 

Note that 

|Z(*)l<Hooy rh dr |pg ,_ p , g| (0-3) 

oo \ e 2iwup V ^ + \ e ^QV(P)\ + f 00 \QV(Q)\ 



\PQ'-P'Q\ J rh \PQ'-P'Q\ 

We will show below that each of these four integrals scales as 1/\oj\ in the large uj limit, 
provided Imw > 0: 



iPQ'-P'Ql ~ \u\ J rh ■ |PQ'-P'Q| ^ |w| 
00 dr/ |e 2 ^QP(P)| < Cqp ^ \QV(Q)\ < Cqq 



|PQ'-PQ| " |w| '7 rh |PQ'-P'Q|- |w| ' 
The inequality (|C.3| ) reduces to 

< (K 2 C PP + K(C PQ + C QP ) + C QQ ) . (C.4) 

Thus if K is of order one compared to u, then Z(s) G D. 
Next we evaluate the l.h.s. of equation (|( 



|J(s) -X(s )| < ||s - solloo ( \\s + solloo / dr' 



r , , \e 2i ^{PV(Q) + QV(P))\ \ 

Jr h \PQ'~P'Q\ J 



PV(P) 



PQ> - P'Q 



(C.5) 



< ( 2K Cpp + c PQ + c QP ) lls _ sq||oo 

M 



is a contraction mapping with c = 2KCpp + Cpq + Cqp. 
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poo 


e Aiuu PV(P) 


poo 

dr = 


e 4iu)U PV(P) 




PQ 1 - P'Q 


Jo 


PQ' - P'Q 



A Lemma 

To demonstrate that 

Fdu < ^ , (C.7) 

we first perform a change of variables so that the integral is over u. Then we split the 
integral. The Hankel functions H>^ +1 y 2 (x) have the property that \H^ +1 y 2 (x)\ = 0(x~ 1 / 2 ) 

as x —7- oo, \H^ +l y 2 { x )\ = 0(x _ ( n+1 )/ 2 ) as x —7- 0. The quantity y(u) has the property that 
|y(tt)l = 0(u _1 ) as u — > oo. We assume that y(u) = 0(u 2A_1 ) as u — > for some A. We 
split the integral into three parts according to whether u is greater or smaller than some fixed 
value uq, and according to whether \oj\u (the absolute value of the argument of the Hankel 
function) is greater or smaller than one. So the three regions are (uo,oo), (|w| _1 ,iio), and 

(o,^- 1 ). 

On (uo,oo) we get a bound of O((\oj\u)~ 1 )0(u~ 1 )du = Odu^ 1 ). On , Uq) we 

get a bound of O^^u^O^-^du, which is O(\oj\- X ) if 2A > 1. On (0, \uj\~ 1 ) we 

get a bound of f^ 1 '' 0((\u\u)- n - 1 )0(u 2A - 1 )du, which is 0(\u\- 2A ) if 2A > n + 1. So, for 
n > 0, 2 A > n+ 1, the integral converges and is 0(|w| _1 ) for large uj. 

For the current-current correlator, we have n = d—3, so we have convergence if 2A > d—2. 
For the stress tensor correlator, n = d — 1, so we have convergence if 2 A > d. 

The result for Cpq, Cqp, and Cqq and the other three integrals follow analogously. So, 
assuming A is sufficiently large, si—)- J £(s) is a contraction mapping with contraction factor 
0(|w| _1 ). We can use this mapping to find an asymptotic expansion in powers of for 



the solution to ( 4.18 ). 
C.2 Fermions 

We show that l(s) Q5.18 ) is a contraction-like mapping. (We assume Imw > 0.) 

^ =1(1(1(0))) -1(1(0)), (C.8) 

^U) = i°0+ 1 \o) -l oj (0) (C.9) 
where T j (0) = 1(1(...1(0))). Equivalently, X^'(0) = V (0) + + ••• + ■ These con- 

3 

tributions are the corrections, which extend the solution ij)^ to nonzero momentum k and 
introduce a non-vanishing scalar <I>. If ^ conver g es > then it converges to a solution of 



(PD 



We would like to show that iJ)W decreases exponentially with j, which would imply that 
the sum converges. First, we observe from the definitions of 1 and J that l(ipi) — l(ip2) = 
J(ijj\) — J(^2) = J(i>\ — { 4>2)- We apply this identity to the right hand side of ( |C.9|) to obtain 



^0') = J(l°i(0) -1^-^(0)) = J^' 1 )). (C.IO) 
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It turns out that is not necessarily smaller than ip. However, we can still show that 

<J{<J(ip)) is smaller than ip by a factor of 0(u~ min ( 1 ' A )) ) where A is the scaling dimension of 
<£. Given our assumption A > 1, we have that min(l, A) = 1. We write 



du' ! du" 



u u 



du J(u, u') ip , 



(C.ll) 



(C.12) 

where we have changed the order of integration and defined a new function J(u,u'). 
Another Lemma 

We now examine the properties of J{u, u') defined in equation ( |C.12[ ) in order to prove the 
fact that the corrections converge, i.e. that 

11^*) ii < c^r 1 ii v> (i ~ 2) ii 



(C.13) 



is satisfied for some C, where we define 
Consider 



sup e 



J(u,u) 



J-f2^yt.(u)(u+u')+(v(u)+v(u')) 



du » e -2H^u»+v(u»)) ^_^E ki * + i$ ) y ^£\ \ u/i (C . 14) 



Once again we have a rapidly oscillating integral. We can bound the integral using integration 
by parts. Let / be the integrand in ( |C.14j ) . 



e iuiu f(u")du" 



1 

M 



ie wu " f{u")\ u u + i I e^ u f'(u")du" 



< |u;r 1 ^2sup|/(n")|+ / \f'(u")\du' 
= OiU- 1 ). 



(C.15) 



since / and /' decrease exponentially as u becomes large. The quantity on the first line of 
( |C.14 ) decreases exponentially as v! becomes large. Then 



J(u, u')tp{u')du' 



We have thus shown that ^ ip^ converges to a solution of the equation of motion 

3 



0(|w^ 1 e- u, )V'(n / )dn , 



(C.16) 
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D. On the Absence of Branch Cuts 



In field theory, retarded Green's functions typically have branch cuts in addition to poles in 
the lower half of the complex frequency plane. Because of our reliance on classical gravity and 
backgrounds with non-extremal black hole horizons, the expectation is that the field theories 
are in a large N strong coupling limit and at nonzero temperature. In this case, the common 
lore is that the Green's functions will have only poles p4| , [45| ], Indeed, the exact results (at 
T 7^ 0) we find in Appendix A.l| lack branch cuts. 22 

On the gravity side, Gn(oj,k) is free of branch cuts in the upper half of the complex 
frequency plane essentially because the differential equation ( |2.7| ) is holomorphic in oj. In the 
next section, we argue that given an assumption about the singular points of the differential 
equation (2.7), the singularities of Gu(u,k), possibly away from a set of discrete points on 
the negative imaginary axis, are entirely determined by the quasinormal mode solutions we 
studied above. The assumption is the requirement that r = and r = oo be regular singular 
points of (|2.7|), and that there be no other "nearby" regular singular points. More specifically, 
we want the Frobenius power series solutions at r = and r = oo to have an overlapping 
region of validity along the real line < r < oo. 23 



We repeat our analysis for fermionic Green's functions in Appendix D.2 , In Appendix [ 
we outline a method that will allow for more general types of singularities in the differential 
equation. This second method does not rule out branch cuts in the lower half plane, and may 
be useful for studying gravity systems away from the large N and strong coupling limit. 



D.l Bosons 

Given our assumption, there exists a Frobenius series solution to ( [2.7D at r = oo which is 
holomorphic in oj on the interval Th < t < oo. The series solution at r = oo will generically 
take the form: 

oo oo 

A(r) = ar^ 1 - 1 ^ aj r~ j + 6r~ A ^ b jr ~ j (D.l) 

3=0 3=0 

where —A and A — n — 1 are called indicial exponents and satisfy the quadratic relation 
A(A — n — 1) = m 2 L 2 . (We are allowing Y to include a mass term.) The coefficients a and 



b are arbitrary and independent of r, and we choose an = bo = 1. The fact that (D.l) is 
holomorphic in oj is perhaps obvious: The indicial exponents of the series do not depend on 
oj, and therefore a,j and bj will be polynomials in oj. Note that if 2A — n — 1 is an integer, 
the expansion may be modified to include logarithms. But the logarithms are functions of r, 
not oj. 



22 In regulating the large frequency behavior of these exact Green's functions, we often add logarithmic terms. 
However, we are free to choose the branch cut of the regulator to lie in the lower half plane. 

23 This assumption can probably be weakened and the following argument made to work provided the 
"nearby" singular points are regular and their indicial exponents independent of uj. 
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The existence of a holomorphic series solution at r = is less obvious as the indicial 
exponents depend on to. We have: 



OO 



A±{r) = ( r _ rft )±W^ ^ c±J .( r _ rh )i . (D.2) 

i=o 

We assume c±o = 1- The coefficients of the ingoing Frobenius series ^4_(r) will involve poles 
at regular intervals 2iw/Ff l = n where n is a positive integer. The poles come from the fact 
that for these values of u, the two horizon power series A±(r) solutions overlap and A_ 
should be modified to include ln(r — r^) dependence. There is a similar problem with u = 0. 
Nonetheless, we can conclude that A- (r) is holomorphic away from these particular values of 
cu and that moreover A- (r) is holomorphic in the upper half of the complex frequency plane. 



To determine a and b in (Chi), and hence the Green's function, we can match the ingoing 
Frobenius power series at the horizon ( p.2[) (and its derivative) to the one at the boundary 
( p.l[ ) at some intermediate point < r < oo. The matching involves solving a linear system 
of two equations and two unknowns, a and b. Because there must exist a globally well defined 
solution, the determinant of the system cannot vanish, and the resulting a and b must be 
holomorphic in u, away from the special values mentioned in the previous paragraph. 24 

The Green's function can be computed from the near boundary limit of A' /A where we 
know a and b are holomorphic in the upper half plane and at worst meromorphic in the lower 
half plane. Hence, any branch cut or singular behavior of the Green's function in the upper 
half plane must come from zeroes of a, which we studied and largely ruled out, except at 
ijj = 0, in the previous subsections. 

If we relax any of the three initial assumptions of large N, strong coupling, or T / 0, 
the argument above generally ceases to work. So, for example there is a possible source of 
confusion regarding the T = limit where we know that the Green's functions do often have 
branch cuts. From the gravity side, the regular singular point at r = can become an ir- 
regular singular point at r = in the limit T — > 0. The asymptotic series expansions around 
irregular singular points typically have a zero radius of convergence, and the matching argu- 
ment above fails. Similarly, 1/N corrections and finite coupling are both known individually 
to introduce branch cuts [f44| , |45[1 in the lower half of the frequency plane. The technique 
outlined in Appendix p.3| may be useful for these cases although it does not work for T = 0. 

D.2 Fermions 



Similar to the discussion of branch cuts in the bosonic case of Appendix D.l, we argue that 
there are no branch cuts in our fermionic Green's functions away from T = 0. Instead of 
working with a second order differential equation, we find it simpler to work directly with the 
first order system (5.7). 



24 The determinant is the Wronskian for the two power series solutions in (D.l) for which Wf = (2 A — n — 
l)r~" / F{r). Superficially then we appear to have missed 2A = n + 1. In this case, an extra logarithm in one 
of the series expansion guarantees that Wf does not vanish. 
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Comparing with the bosonic case of Section D.l , we make the same assumptions about 
the regular singular points of the differential equation, and we attempt to match a power series 
solution at the horizon to a power series solution at the boundary. The fermion equation of 
motion ( |5.7|) may be written 

< = tya (D.3) 

where X is a 2 x 2 matrix encoding the coefficients from ( |5.7| ). In computing these power 
series, diagonalizing the matrix X at the regular singular points may introduce square root 
branch cuts. In fact, the boundary and horizon power series expansions are free of such cuts. 
For the massless fermion, the boundary is actually not singular, and no diagonalization is 
required. Near the horizon, we know 



(r - r h )F, 




(D.4) 



The eigenvalues and eigenvectors of lim r _;. rh (r — r^)T have no square root branch cuts. The 
Frobenius power series for the infalling solution will have at most logarithms and only for 
the same special values of oj along the negative frequency axis that we found in the bosonic 
case. Thus the horizon and boundary expansion of ifi are holomorphic in uj everywhere except 
for possibly a discrete set of uj along the negative frequency axis. We can conclude the only 
singularities in the fermionic Green's function (away from these discrete values) come from 
quasinormal modes. It should be possible to relax the assumptions about the regular singular 



points of (5/7) along the lines of Appendix D.3 



D.3 An Extended Class of Models 

We outline a refined technique for demonstrating the absence of branch cuts in the upper 
half of the complex plane; we consider only Imw > and uj ^ 0. In proving the absence of 
branch cuts before, recall that we restricted our differential equations to have a very limited 



class of singular behavior in Section D.l and p.2j . We believe this restriction is consistent 
with field theories at T ^ 0, large N and strong coupling. The technique here is valid for an 
extended class of models where the Green's functions come from solving differential equations 
with worse types of singularities. Such singularities may occur in moving away from the large 
N and strong coupling limit in field theory. However, this particular technique will fail at 
T = 0. 

We make the following assumptions: 

• p and Y are integrable on the interval (r#,oo), and F is positive. 

• Near the horizon, F ~ r — r# and the integral of Y converges. 

• Near the boundary, F ~ r and Y = 0(r- 2 ). 

Let s = ~ tFl A — , then the differential equation for s is 

is? -4— + + YrH = 0- (D.5) 

P r n F K ' 
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Suppose we write s as a Taylor series about some ujq in the upper half plane: 



s(r,oj) = J2sk(r)(io-uo) k - (D.6) 



k=0 

If the Taylor series has a nonzero radius of convergence, then s is analytic in a neighborhood 
of uq. We can write an infinite sequence of differential equations: 

71 

*4 = E - (^ow 2 + 25 kl oj + 5 fc2 )^ - Yr n 5 k0 (D.7) 



From our analysis in Section 4.3 and Appendix C.l, we know how so behaves. Note that, for 
k > 0, all of the equations are linear in s k . Since the imaginary part of s °pJjP is non-positive, 
there exists a solution for each s k that remains finite at the horizon. We see that the s k 
obey a recurrence relation similar to that of the Catalan numbers. It is not hard to show 
using generating functions that the s k grow at most exponentially, and thus the radius of 
convergence is finite. 

In the limit r — > oo, we have to worry about the s k diverging. Generically, the terms 
so,si,S2 are 0(r™ -1 ). Because of the Fr n in the denominator of ( D.7 ), subsequent terms 



diverge successively more slowly, and for sufficiently large k there is no divergence. So we 
may renormalize s by subtracting a function that is polynomial in u). 

Since J and f -S&n are bounded as r — > oo and all but finitely many s k are bounded 
as r — > oo, our argument that the s k grow at most exponentially works in the limit r — > oo 
as well. 
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